Optimal control, geometry, and quantum computing 
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I. INTRODUCTION 

Quantum computers have caused great interest due 
to their potential use in efSciently solving problems 
considered intractable on conventional classical comput- 
ers Despite this interest, there is as yet no general 
framework for constructing efficient quantum algorithms, 
nor for proving limitations on the power of quantum com- 
puters. 

Recent work 0, has proposed a geometric approach 
to quantum computation, based on the observation that 
finding quantum circuits of the minimal size required 
to perform some desired computation is equivalent to a 
problem in Riemannian geometry. More precisely, the 
size of the minimum quantum circuit synthesizing a uni- 
tary U is, up to polynomial factors and some technical 
caveats (see Section HVI for precise statements), equal to 
the distance d{I, U) between the identity operation / and 
[/, according to some Riemannian metric. This equiva- 
lence means that problems in quantum computation can 
be recast in terms of equivalent problems in Riemannian 
geometry. 

The results of 0, establish an equivalence between 
the number of gates needed to synthesize U and the 
minimal distance according to some specific Riemannian 
metric. However, inspection of the proof in 3, 4J shows 
that many of the properties used in the proof are rather 
generic, and there are certainly other Riemannian met- 
rics with the same property. One may therefore ask what 
is the most general class of Riemannian metrics that can 
be connected to gate complexity. Even more generally, 
the problem of finding minimal geodesies in Riemannian 
geometry may be viewed as an instance of the problem of 
optimizing some cost function in the framework of non- 
linear optimal control (see, e.g. Q), and it is interesting 
to ask whether it is possible to make any general connec- 
tions between optimal control and gate complexity. 

The purpose of the present paper is to identify a large 
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family of optimal control problems whose optimal cost 
is equivalent to the minimal gate complexity of the de- 
sired unitary operation. As special cases of our results 
we obtain the geometric results of 's', , but also iden- 
tify many other classes of optimization problems which 
can be connected to gate complexity, including problems 
from time-optimal control, and from Riemannian, sub- 
riemannian, and Finslerian geometry. Of course, in some 
(though not all) of these examples more straightforward 
techniques may be used to relate the optimal cost to 
quantum gate complexity. The benefit of the analysis 
in the present paper is that it provides a unified and 
generalized framework for deriving connections between 
quantum gate complexity and optimal control. 

By identifying this large family of optimal control 
problems we identify the essential features of the geomet- 
ric problem in [j, |j| that are responsible for the equiva- 
lence to quantum computation. We also widen the class 
of problems in optimal control which may be analysed 
in order to obtain insight into quantum computation. 
A considerable body of work has been done on optimal 
control in quantum physics (see references later in the 
paper), and we hope that the close connection between 
optimal quantum control and quantum gate complexity 
will stimulate further work on optimal quantum control. 

The structure of the paper is as follows. Section ^de- 
scribes background material on quantum computing and 
optimal control theory that is useful later in the paper. 
Section [lIII proves a general theorem relating the optimal 
cost for a control problem to quantum gate complexity. 
In Section IIVI we illustrate this theorem through a se- 
ries of applications to example problems, including time- 
optimal control, and problems from Riemannian, subrie- 
mannian and Finsler geometry. Section IVI concludes . 



II. BACKGROUND 

In this section we introduce some background material 
on quantum computation (Subsection III A|l and optimal 
control (Subsection IIIB|) that will be useful later in the 
paper. 
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A. Quantum computation and gate complexity 

We assume the reader is familiar with basic notions 
of quantum circuits (e.g., Chapter 4 of 2]). Suppose 
U is an n-qubit unitary operation. We define the ex- 
act gate complexity G{U) to be the minimal number of 
one- and two-qubit quantum gates required to synthe- 
size U exactly, with no ancilla qubits allowed to assist 
in the preparation of U. We define the approximate gate 
complexity G(U, e) to be the minimal number of gates 
required to synthesize some n-qubit unitary operation V 
satisfying \\U — V\\ < e, where || • || is the usual matrix 
norm. Once again, no ancilla qubits are allowed to assist 
in the synthesis. Note that in (3, |^] the notation ■m{U) 
was used for the gate complexity. 

Our results connect problems in optimal control to the 
values of G{U) and G{U, e). The typical object of interest 
in optimal control is the optimal cost C{U) associated to 
a unitary, U, according to a cost function which is defined 
precisely below. Our goal is to identify control problems 
such that G{U) provides good lower bounds on the exact 
gate complexity G(U), and good upper bounds on the 
approximate gate complexity G{U,e). As a result, up 
to polynomial factors the exact synthesis of U without 
ancilla must take at least G{U) quantum gates, and U 
can be synthesized to accuracy e using at most C(U) 
quantum gates. 

One might naturally ask if it is possible to extend 
these results to prove a similar lower bound involv- 
ing approximate computation, or an upper bound in- 
volving exact computation. Parameter counting can 
be used to show that a bound of the form G{U) < 
poly{C{U),n) is not possible. Whether a bound of the 
form poly(C([/), n, 1/e) < G{U, e) is possible remains an 
open problem. Fortunately, lower bounds for exact com- 
putation and upper bounds for approximate computation 
remain of great interest. 



B. Optimal control on SU{2") 

We now sketch the basic ideas of optimal control the- 
ory, following the standard approach (e.g., |5||), but omit- 
ting mathematical details regarding smoothness and reg- 
ularity conditions, as these are not important for our 
purposes. 

Let Hi , . . . , Hrn be a set of linearly independent ma- 
trices in the Lie algebra su(2") of traceless n-qubit 
Hermitian matrices^. Our control system is based on 



^ Note that physicists' and mathematicians' definitions of Lie al- 
gebras differ by a factor of i, and so our definition of su{2") 
is consistent with the usual mathematical definition in terms of 
traceless sA;ei«-Hermitian matrics. 



Schrodinger's equation: 

— = -^H{t)U{t)■ H{t)^Y.^,{t)H,, (1) 

where h{t) = {hi{t), . . . , hm{t)) is known as the control 
function, and we impose the initial condition f/(0) = /. 
Defining the notation = X]j=i^j^jj that 
H{t) = Hh{^t)- We refer to H{t) as the control Hamil- 
tonian corresponding to the control function h{t). Note 
that to any control Hamiltonian H[t) defined on an in- 
terval [0, T] there exists a unique solution U{t) to Equa- 
tion defined on the same interval. 

In general, the control function h{t) is not allowed to 
take arbitrary values, but is constrained to lie in an al- 
lowed control region A C i?™. We denote the correspond- 
ing set of allowed control Hamiltonians by H-a- 

To complete the specification of the control problem we 
must also specify a cost function, which is a real-valued 
function c : A ^ R ow the allowed control region. Equiv- 
alently, it may be regarded as a function c : Ti.A — > R 
on allowed control Hamiltonians, and it is this view- 
point we shall take most often. The cost function allows 
us to assign a cost to a control Hamiltonian H{t) de- 
fined on an interval [0,T] by C{H{t)) = dtc{H{t)). 
This allows us to define the cost of a unitary U by 
C{U) = inix.H{t) G{H{t)), where we take the infimum 
over all intervals [0, T] , and over all control functions H{t) 
such that H{t) e Ha for aU times t, and U{T) ^ U. Note 
that in general there is no reason why this infimum should 
exist, as there may be no allowed control Hamiltonian 
H{t) which can be used to synthesize the desired unitary 
U. However, if we assume that the Lie algebra generated 
by Hi, ... , Hm is the full Lie algebra sm(2"), and that the 
allowed control region TCa is not trivial, we can ensure 
that such a control function exists, and so the infimum 
is defined |^. This condition is known as the condition 
that the control system be bracket-generating. Provided 
reasonable continuity assumptions are made about the 
cost function c(-) it can also be shown that the infimum 
is achieved for some control function H{t). 

The allowed control region T-Ca and the cost function 
c(-) jointly specify the control problems we shall be in- 
terested in. Such control problems are known as right- 
invariant control problems on the Lie group S'C/(2"), and 
we shall denote them using the notation {Ha, c). 

III. BOUNDS RELATING OPTIMAL CONTROL 
AND QUANTUM GATE COMPLEXITY 

In this section we develop some general relationships 
between the cost function C{U) of a right-invariant con- 
trol system {Ha, c) on SU (2") and the exact and approx- 
imate gate complexities, G{U) and G(U,e). Our results 
generalize and extend the ideas in 

Splittings: The key tool we use to relate the cost 
C{U) to the gate complexities G{U) and G{U,e) is an 
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object we refer to as a splitting. We define splittings in 
two steps. First, we identify a special set Hp C Ha 
of preferred Hamiltonians, which we shall assume are 
bracket generating. Second, we identify a projection 
map P : Ha Hp which takes any allowed Hamil- 
tonian H and projects it onto a preferred Hamiltonian 
Hp = P{H). Note that this can be an arbitrary func- 
tion, and need not be a projection in the linear algebraic 
sense. We call the pair [Hp , P) a splitting for the control 
problem [Ha, c). 

The bounds relating the control cost C{U) to gate com- 
plexity will depend on the particular splitting we choose. 
For examples of "good" choices of splitting (i.e., choices 
resulting in fairly tight bounds between control cost and 
gate complexity) see the later examples. For now we sup- 
pose that the choice of splitting has been fixed, and will 
show how it can be used to relate the control cost to gate 
complexity. 

Our construction is rather abstract, and many readers 
may prefer to first read the statement of Theorem^ and 
then to read section HVI where that theorem is applied 
to several example control problems. 

Relationship between C{U) and G{U): To ex- 
press this relationship we need to define two quanti- 
ties associated to the splitting. The first quantity is 
the maximal cost of applying any preferred Hamiltonian, 
cp = sup^g^^ c(iJ). Note that we use the subscript P 
as a mnemonic to indicate that cp is a cost associated to 
the set of preferred Hamiltonians. The second quantity is 
the maximal time Tp required to exactly generate an ar- 
bitrary one- or two-qubit unitary operation by applying 
time-dependent preferred Hamiltonians. 

Observe that we can synthesize any one- or two-qubit 
quantum gate for a cost at most cpTp. Since U can be 
synthesized exactly using G{U) one- and two-qubit gates, 
we deduce the desired bound relating C{U) and G{U): 

G{U) < cpTpG{U). (2) 

Note that the value of G{U) depends only on the control 
system, {Ha,c), not on the choice of sphtting, {Hp,P). 
Thus, different choices of splitting can give rise to dif- 
ferent bounds, and it is necessary to choose the splitting 
in an intelligent way to get the best possible bound. In 
particular, one should choose the splitting to minimize 
the product cpTp. 

Relationship between C{U) and G{U, e): This rela- 
tionship is rather more complex than that between C(U) 
and G(U), and is expressed in terms of four quantities 
associated to the splitting. The first quantity is a ra- 
tio defined by^ R = maxnenA \\H - Hp\\/c{H). The 
second quantity is the maximum matrix norm Np = 
max^fg-Hp of any preferred Hamiltonian. 



^ Note that here and elsewhere we write max and min rather than 
sup and inf. Our proofs are easily modified for the case when 
(for example) the maximum is not defined, but this does make 
the discussion less transparent, and so we have avoided it. 



The third quantity requires a more complex expla- 
nation. Suppose A > and (5 > 0. We define a A- 
averaged Hamiltonian to be a Hamiltonian H which can 
be written in the form H = dt H{t) for some Hamil- 
tonian control function which remains in the preferred 
set, H{t) & Hp. We define the A-averaged unitaries to 
be the set of unitary operations which can be written in 
the form exp(— ii?) for some A-averaged Hamiltonian H. 
We define g{ A,S) to be the maximum number of one- and 
two-qubit gates required to approximate an arbitrary A- 
averaged unitary to an accuracy better than d in matrix 
norm. 

The fourth quantity is the minimal cost associated to 
any allowed Hamiltonian, ca = niinj^gT^^ c{H). This 
quantity arises in our proof as a way of getting a bound 
on the time T associated to the optimal Hamiltonian 
control H{t). The argument is to observe that C{U) = 

dtc{H{t)) > TcA, and so T < G{U)/ca. 
With these quantities defined, we can relate C{U) and 
G{U,e). The first step is to take the Hamiltonian con- 
trol H{t) which achieves the optimal control cost C{U), 
and to form the corresponding projected Hamiltonian 
Hp{t) = P{H{t)). We suppose Hp{t) generates a uni- 
tary Up, and aim to show that C/p is a pretty good ap- 
proximation to U . As in the proof of Lemma 1 in the 
supporting online materials for Q , we can apply the tri- 
angle inequality repeatedly to obtain: 

\\U-Up\\ < [ dt\\H{t)-Hp{t)\\. (3) 
Jo 

The definition of the ratio R ensures that \\H — Hp\\ < 
Rc{H) for all H, and thus: 

j dt\\H{t) - Hp{t)\\ <R [ dtc{H) = RC{U). (4) 
Jo Jo 

Putting these inequalities together we obtain ||[/ — t/pjl < 
RC{U). Intuitively, provided the control problem and 
splitting are such that R is much smaller than l/C(/7), 
we ensure that U and Up will be quite close. 

In the next step of the proof we discretize the evolu- 
tion according to Hp{t), and show that it can be approx- 
imated by a suitable sequence of A-averaged Hamiltoni- 
ans. The key to doing this is the following lemma, which 
appeared as Lemma 2 in 0. We have made some minor 
notational changes to the statement of the lemma, but 
the essential content of the lemma, and the proof, which 
is an easy application of the Dyson operator expansion, 
is unchanged. 

Lemma 1 Let V be an n-qubit unitary generated by ap- 
plying a time- dependent Hamiltonian Hp{t) £ Hp over a 
time interval [s, s -|- A]. Then defining the corresponding 
A-averaged Hamiltonian H = J^~^^ dtH{t) we have: 

\\V - e-'"^\\ < 2(e^^^ - 1 - iVpA) = 0(7V|,A2), 

(5) 
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where Np is the maximum matrix norm of any preferred 
Hamiltonian, as defined earlier. 

To apply this lemma, we divide the time interval [0, T] 
up into a large number N of time intervals each of length 
A = T/N. Let Up be the unitary operation generated 
by Hp{t) over the jVci time interval. Let C/^^ (the unitary 
corresponding to the mean Hamiltonian) be the unitary 
operation generated by the A-averaged Hamiltonian over 
the corresponding time interval. Then the lemma im- 
plies that ||[/^ - Ul^W < 0{NlI\^). By assumption, 
we can then synthesize a unitary operation Uj^ using 
at most g(A, 6) one- and two-qubit gates, and satisfy- 
ing II — C^ill < ^- We define Ua (the actual unitary 
to be synthesized by our gate sequence) to be the result 
of applying the unitaries C/^ in sequence. Note that Ua 
can be generated using Ng{A, S) = Tg{A, S)/A one- and 
two-qubit quantum gates. 

Repeated application of the triangle inequality, substi- 
tution of the inequalities obtained above, and using the 
fact that N ^T/A, yields: 

\\U-Ua\\ (6) 

< II {/- C/p 11 + II f/p-C/A II (7) 

N 

< RC{U) + J2\\Uj,-U\\\ (8) 

N 

< RC{U) + J2[\\Uj,-Ui,\\ + \\Ui,-Ui\\) (9) 



< RC{U) + 0{NI,TA) + —S. 



(10) 



Substituting the bound on T obtained earlier, T < 
C{U)/ca, we deduce that we can synthesize an opera- 
tion Ua satisfying 



r-t/.ii<«c(t/)+o|^4W2^)+S2^ 

CA J caA 



(11) 



using C{U)g{A,6)/cAA one- and two-qubit gates. 
Summing up, we have the following theorem: 



Theorem 1 Consider a control problem {Ti.A,c) and a 
splitting (Tip^P) for that problem. Then we have: 

(1) Let cp = niaxHeHp c{H) be the maximal cost 
of any preferred Hamiltonian, and suppose Tp is the 
maximal time required to generate an arbitrary one- or 
two-qubit unitary operation using preferred Hamiltoni- 
ans. Then: 



C{U) < cpTpG{U). 



(12) 



(2) Let R = maiCHeHA \\H - Hp\\/c{H), Np ^ 
maxHeT-ip \\H\\, ca = 'nrniHeHA ^{H)- Suppose that if H 
is a A-average of Hamiltonians in Tip, i.e., can be writ- 



control function H{t) which remains in the preferred set, 
then the corresponding unitary exp(— iff) can be simu- 
lated to an accuracy S using a number of gates g{A,6). 
Then we can synthesize an operation Ua satisfying 



\\U-Ua\\ <RC{U) + 



NlC{U)A \ C{U)5 

CA ) CaA 



(13) 



using C{U)g{A^5)/cAA one- and two-qubit gates. 



ten in the form H — dtH{t) for some Hamiltonian 



We stress that this theorem does not necessarily give 
tight connections between optimal costs and gate com- 
plexity. Finding such connections depends on making 
an appropriate choice of the cost function, and of the 
splitting. However, the examples in the next section will 
show that such choices can be made for a wide variety of 
interesting cost functions. 



IV. EXAMPLES 

We will now describe a sequence of examples illustrat- 
ing Theorem ^ These examples are not exhaustive, but 
illustrate the wide range of situations in which Theo- 
rem ^ can be used to relate problems of optimal control 
and quantum gate complexity. 

Note that in each of the examples described in the 
present section, we are imagining that there is a fam- 
ily [/ = t/„ of unitary operations, one for each value of 
n, acting on n qubits. Correspondingly, in each of our 
examples we will describe an entire family of cost func- 
tions and splittings, one for each value of n. Our goal 
is to prove results of the form poly(C(C/), n) < G{U) 
and G{U, e) < poly(C(f7), n, 1/e) for suitable polynomial 
functions. 

Subriemannian metric: Suppose the allowed Hamil- 
tonians Ti.A are of the form H = '^^ha-cr, where the 
sum is restricted to be over Pauli sigma matrices con- 
taining only one- and two-qubit terms, and we require 
that J2a — ^- We define the cost function by c{H) = 
VScr f^a- so for allowed Hamiltonians we have c{H) — 1. 
This cost function C{U) is an example of the distance 
associated to a subriemannian metric 1^, and the prob- 
lem of finding C{U) is that of finding the minimal length 
geodesies on a subriemannian manifold. We choose the 
splitting to be trivial, with Tip = Ha and P{H) = H . 

With this control problem (?iyi,c) and splitting 
(TLp,P), we may apply part (1) of Theorem 1. In that 
notation, it follows immediately from the definitions that 
Cp = 1 and Tp is a constant of order one, independent 
of the number of qubits, n. Thus C{U) < TpG{U), and 
so, up to a constant factor, the subriemannian distance 
C{U) provides a lower bound on the exact gate complex- 
ity G{U). 

To apply part (2) of Theorem ^ note that we have 
R ~ Q and ca — again directly from the defini- 
tions. It follows from elementary norm inequalities that 
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Np = 0{n)^. To understand the behaviour of g{A,S), 
suppose that H = J^dtH{t) is a A-averaeed Hamilto- 
nian over Hamiltonians in Tip. Lemma 3 in JJ imphes 
that exp{—iH) can be simulated to an accuracy of order 
0{n^A^) using 0{n^/A) gates. Thus g{A,0{n'^A^)) < 
0{r?IA). We deduce that we can synthesize an opera- 
tion Ua satisfying 

\\U ~Ua\\< 0{C{U)n'A) + 0(C{U)n^A^) (14) 

using 0{C{U)n? / A"^) gates. It follows that by choosing 
A appropriately, we can synthesize a good approximation 
to U using a number of gates that scales in a fashion 
comparable to C{U). To see this, let A — e/n^C(C/). 
Then we see that we can synthesize an operation Ua 
satisfying \\U ~Ua\\ < 0(e) using 0{C{JJfn^ / e^) gates. 
It follows that: 

G{U, e) < 0{C{Ufn^le^), (15) 

which is the required result — G{U, e) scales as no more 
than a polynomial in C{U), n and 1/e. 

Time-optimal control: If c{H) = 1, then C{U) 
is the minimal time taken to generate U using control 
Hamiltonians in the allowed control region, Ti.A- This 
is known as the time-optimal control problem. A com- 
mon variant of the time-optimal control problem is to 
constrain the set of allowed controls so that hi{t) = 1, 
i.e., so that the Hamiltonian Hi is always being applied. 
This is known as the time-optimal control problem with 
drift, and Hi is known as the drift Hamiltonian. The 
time-optimal control problem in quantum physics has re- 
ceived considerable attention; see, e.g., lol ITol lTllIT^ 
for recent work, and further references. Of particular 
interest in this context is work such as Q, which stud- 
ies the time complexity of various quantum computing 
primitives, such as the quantum Fourier transform, and 
applies powerful tools from optimal control theory such 
as the Pontryagin maximum principle |l3j| (see, e.g., 0) 
to obtain time-optimal implementations of these primi- 
tives. 

The time-optimal control problem with drift takes a 
particularly simple and appealing form in the case where 
there are only two terms in the control Hamiltonian, 
i.e., H — Hi + h{t)H2, and it is this case we shall fo- 
cus on; analogous results can also be proved for other 
time-optimal control problems using essentially the same 
ideas. We will assume that the control region is such that 
the allowed range of values for h{t) is \h{t)\ < I. A priori 
it is not obvious that it is possible to find examples of 



Hamiltonians Hi and H2 which are bracket-generating. 
However, it follows from results of p^lTH (c.f. ^3|) that 
if we choose Hi and H2 at random, then with probabil- 
ity one they will be bracket-generating. Of course, this 
does not mean that they are universal for quantum com- 
putation in the usual sense. It may take such a Hi and 
H2 exponential time to generate standard quantum gates 
such as the controUed-NOT, or even single-qubit unitaries. 
Conversely, it may not be possible to efficiently simulate 
Hi and H2 in the standard quantum gate model of com- 
putation. 

We will now provide examples of families of Hamilto- 
nians Hi and H2 such that the time-optimal control cost 
scales as a polynomial in the quantum gate complexity. 
The key to this is the following theorem, which is of in- 
dependent interest: 

Theorem 2 There is a family of n-qubit Hamiltonians 
Hi and H2 such that: (1) any one- or two-qubit unitary 
gate can be synthesized exactly in a time bounded above 
by a value that scales as a polynomial in n; and (2) using 
one- and two-qubit gates we can simulate any unitary 
of the form exp(— iA(_ffi -f aH2)) (with \a\ < \) to an 
accuracy S using g{A,d) = 0{p{n)A'^ /6) one- and two- 
qubit gates, for some polynomial p{n). 

Proof (outline): We choose Hi to be a Hamiltonian 
acting on the first two qubits in a manner specified more 
precisely below. We choose H2 so that exp(— 1^2) per- 
mutes qubits 2 through n by a cyclic displacement, i.e., 
the state of qubit 2 becomes the state of qubit 3, the 
state of qubit 3 becomes the state of qubit 4, and so on, 
with the state of qubit n becoming the state of qubit 2. 

With these choices, conclusion (2) follows from stan- 
dard quantum simulation techniques for simulating a sum 
of Hamiltonians, and the observation that the Hamilto- 
nians Hi and H2 can both be efficiently simulated (the 
latter using the quantum Fourier transform P,I3)- 

Conclusion (1) requires a little more effort. In partic- 
ular, note that using Hi and H2 we can simulate the 
Hamiltonian cxp{~iH2)Hiexp{iH2) = Hi, where the 
tilde denotes that Hi is the same Hamiltonian as Hi , but 
now acts on qubits 1 and 3. It can now be verified nu- 
merically or by hand that for many choices of two-qubit 
Hamiltonian Hi, the Hamiltonians Hi and Hi generate 
the full Lie algebra on qubits one, two, and three''. As a 
result, in constant time we can generate an arbitrary uni- 
tary operation on qubits one, two and three. Conjugating 
repeatedly by exp(— 2^2) we can use this to generate an 



3 ll^^ll < Eal'^-I ^ (V3/2)n^E^ft2 = {v^/2)n, e Hp. 
The second inequality follows from \\v\\i < v^||'Ci|2 where d is 
the dimension of the real vector v, \\v\\i = 1^*1 ^^'^ — 

our case d = 9n(n — l)/2 + 3n, the number of 
one- and two-qubit terms. 



^ Examples of this phenomenon were f oun d numerically by the 
presen t au thor and H. L. Haselgrove lia. C. Hill and Hasel- 
grove [ij have recently constructed rather more elegant ex- 
amples demonstrating essentially the same phenomenon as de- 
scribed in this theorem, but making use of a relatively simple 
(and more physically plausible) two-body Hamiltonian in place 
of H2, which involves complex many-body terms. 



6 



arbitrary unitary on qubits 1 and j, where j is any qubit. 
Standard techniques then suffice to efficiently generate an 
arbitrary unitary on any pair of qubits. QED 

Suppose we consider the time-optimal control problem 
where Hi and H2 have been chosen as in Theorem |5] 
As in the subriemannian case we again choose the trivial 
splitting, Hp = Ha and P{H) = H. Applying part (1) 
of Theorem^ we see that cp = 1 and Tp < q(n), for some 
polynomial q{n). As a result, we have C{U) < q{n)G{U). 

Applying part (2) of Theorem ^ we have i? = 0, 
Np = 0(1), CA = 1 and g{A,5) = 0{p{n)A'^/5), for 
some polynomial p{n). As a result, we conclude that it 
is possible to synthesize a unitary Ua satisfying 



\U-Ua\\ < OiC{U)A)+0 



(16) 



using 0{C{U)p{n)A/6) gates. Setting A = e/C{U) and 
6 — e^/C([/)^, we see that we can synthesize a unitary 
Ua satisfying \\U - Ua\\ < 0{e) using 0{C{U)^p{n)/€) 
gates, and so we conclude that 



G{U,e)<OiC{U)^p{n)/e) 



(17) 



which is the desired polynomial scaling. 

Riemannian metric: We now analyze the metric 
considered in P|, and show how to recover the results 
of Q. This is our first example which makes use of a 
nontrivial splitting. Expanding the control Hamiltonian 
as H = hcr<y, where the sum is over all n-qubit Pauli 
matrices, the cost function of 4] (which is just the norm 
associated to the metric) is defined by: 



\ 



J2hi+p^j2h,, 



(18) 



where the primed sum is over one- and two-qubit Pauli 
terms, and the double primed sum is over three- and 
more-qubit Pauli terms. The parameter p is a penalty 
whose value we set later. Ha is defined to contain all 
those Hamiltonians such that c{H) = 1. For the split- 
ting, we choose the set of preferred Hamiltonians Hp 
so that it contains all Hamiltonians containing just one- 
and two-qubit terms. The projection P takes an arbi- 
trary Hamiltonian, and eliminates all terms in the Pauli 
expansion except the one- and two-qubit terms, i.e., it 

takes X^CT ^cr'^ + Sct ^'^'^ to X^CT ''•fO'- In the language of 
part (1) of Theorem n we have cp = 1 and Tp is a con- 
stant, and so C{U) < TpG{U), i.e., the control cost is a 
lower bound on the gate complexity G(U), to within a 
constant factor. 

Next, we evaluate the quantities defined in part (2) of 
Theorem ^ To evaluate R, observe that the Hamilto- 
nian H achieving the maximum must contain only terms 
which are three- or more-body. Thus: 



2" 



(19) 



where the first inequality follows from the triangle in- 
equality, and the second inequality follows from the 
Cauchy-Schwarz inequality. For the same reasons as 
in the subriemannian case, Np = 0{n), ca = I and 
<?(A,0(n4A3)) < 0{n^/A). 

Applying part (2) of TheoremQlwe deduce that we can 
synthesize an operation Ua satisfying 

on 

\\U-Ua\\< —C{U) + 0{C{U)n'A) + 0{C{U)n^A^) 
V 

(20) 

using 0{G{U)'n? / A^) gates. Again we choose A — 
e/n?C{U). Then we see that we can synthesize an op- 
eration Ua satisfying \\U - Ua\\ < 2"C(C/)/p Cl(e) us- 
ing 0(C(J7)'^ri^/e^ gates. Standard results on univer- 
sality (see, e.g., [13 and references therein) imply that 
C{U) < 0(4") for all unitaries U, and so by choosing 
p = 8"/e we obtain 



G(C/,e) < 0{GiUynye'), 



(21) 



which is the desired polynomial scaling. 

Other control problems: It is not difficult to gen- 
erate many other examples of optimal control problems 
whose cost scales in essentially the same way as the gate 
complexity. An example is the following cost function 
that was conjectured in Q to be equivalent to the gate 
complexity: 



^/iffCT = \hcr\ +py^ \ho 



(22) 



Once again, p is a penalty parameter that we shall choose 
to be large. We define Ha to consist of all Hamiltonians 
such that c{H) = 1. We define the splitting as for the 
Riemannian metric considered above, setting Hp to be 
those Hamiltonians in Ha containing only one- and two- 
qubit terms, and the projection P to remove all three- 
and more-qubit terms from the Pauli expansion. A sim- 
ilar analysis to the Riemannian case allows us to relate 
the cost G{U) to the gate complexity. The only signifi- 
cant difference is in the evaluation of R, where we obtain 
R ^ l/Pi and thus it is possible in this case to choose 
more modest values of p and still achieve a close rela- 
tionship between the scaling of the cost and of the gate 
complexity. 



V. CONCLUSION 

We have proved a general theorem relating quantum 
gate complexity to the optimal control cost for an arbi- 
trary control problem. Application of the theorem de- 
pends on the use of a tool known as a splitting, which 
must be chosen appropriately in order to obtain good 
bounds. We have illustrated this theorem with exam- 
ples showing that quantum gate complexity is essentially 
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equivalent to the optimal control cost for problems in- 
cluding time-optimal control and finding minimal dis- 
tances on certain Riemannian, subriemannian, and Fins- 
lerian manifolds. It is possible to improve the scaling 
in many of these results with a more refined use of the 
Dyson operator expansion 20J and Suzuki- Trotter type 
formulas and it would be interesting to determine 
what the optimal bounds are. It also seems likely that 
the results can be further generalized using tools more 
sophisticated than the notion of a splitting that we have 
introduced. However, the most important direction of fu- 
ture work will be to better understand the optimal cost 



for specific choices of control problem, and what it im- 
plies for quantum gate complexity. 
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